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Abstract We look for singlevalued solutions of the squared modulus M of the 
traveling wave reduction of the complex cubic-quintic Ginzburg-Landau equa- 
tion. Using Clunie's lemma, we first prove that any meromorphic solution M 
is necessarily elliptic or degenerate elliptic. We then give the two canonical de- 
compositions of the new elliptic solution recently obtained by the subequation 
method. 
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1 Introduction. The CGL5 and CGL3 equations 

When a system is governed by an autonomous nonlinear algebraic partial dif- 
ferential equation (PDE), it frequently admits permanent profile structures 
such as fronts, pulses, sinks, etc |28j . and usually these profiles are mathemat- 
ically some singlevalued solution of the traveling wave reduction (x, t) — > x — ct 
of the PDE to an ordinary differential equation (ODE). 
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When the field is a slowly varying complex amplitude A, the simplest 
equation involving time evolution, dispersion, nonlinearity and forcing is the 
one-dimensional complex Ginzburg-Landau equation 

iA t + P A xx + q\A\ 2 A + r\A\ i A-i 1 A = 0, {A, P ,q,r)eC, 7 el. (1) 

We only consider in this class the equations which have the worst singular- 
ity structure, the cubic one (CGL3, r = 0, lm(q/p) ^ 0) and the cubic-quintic 
one (CGL5, Im(r/p) ^ 0). For a summary of results, see the reviews [Tll2"5]. 

CGL5 depends on seven real parameters. Its travelling wave reduction 



/2 



A(x,t) = V / M(£)e l( - w< + ( ^)), £ = x-ct, (c,cj,M,ip) € R, (2) 

// ;2 ,M' c M' c , er , r , r9 
icp" -<p" + up'— - i — — + V + -M + -M 
M 2p M p p p 

■ = 0, (3) 



P 

depends on eight real parameters, denoted e r , ej, d r , di, s r , Sj, g r , gt, 

p p p 

7 + iu c 2 s r ,„ . . . 

9r + Wi = 1 7-{*Si + ls r)- (4) 

p 4 

The fact of taking account of the phase invariance [55] 

<p' = ^j + Re^-, (5) 
2p 

further reduces them to the seven real parameters e r ,ei,d r ,di,g ri gi, csi = 

This third order system ^ can be written either as a real two-component 
rational system in the real variables (M, ib), 

M" M' 2 M' 2 - , , r 

— — — - — — r - Ki— — - ib 2 + e r M 2 + d r M + g. t = 0, 
2M AAjl 2 2M v y (6) 

V 1 ' + i>-^ - Kiip + e t M 2 + diM - g r = 0, 

or, by elimination of ip, as a real third order second degree ODE in M [50] 
I 2njG — G' 2 _ G 

V ~ 2M 2 {e l M 2 + diM — g r ) ~ M 2 ' U 
(G' - 2k[G) 2 - AGM 2 (e t M 2 + d t M - g r f = 0, (8) 

G = -MM" - ^M' 2 - jMM' + e r M A + d r M 3 + gi M 2 . (9) 

The purpose of this work is to show that, for all values of the seven pa- 
rameters, all meromorphic particular solutions M of CGL3 and CGL5 belong 
to class W (like Weierstrass), defined as consisting of elliptic functions and 



Meromorphic solutions of CGL5 



3 



their successive degeneracies, i.e.: elliptic functions (meromorphic doubly pe- 
riodic), rational functions of one exponential e k ^, k € C, rational functions of 
£. The assumption M meromorphic implies the same property for the variables 
M'/M,ip and the logarithmic derivative of the complex amplitude Ae lult , 

B := — log f Ae lujt ~ ics r£/2^ = ML +i ^, ( W ) 
dC V J 2M v y ' 

In case qui = 0, the real system ([5]) displays a parity invariance, 

q = 0: (Af,V,0->(-Af.^0, (11) 
Ki = 0: (M,V,0->(^-^-0- (12) 

This paper presents results complementary to those of [5] . Section [5] recalls 
the singularities of M and ip. In section [31 we prove that, for all values of 
the seven real parameters, any meromorphic solution M of ((5]) is in class 
W. In section [H we recall a method to obtain all the elliptic or degenerate 
elliptic solutions M, and present the first order ODE for M(£) characterizing 
the recently obtained elliptic solution [5] . Finally, section [5] is devoted to the 
construction of canonical expressions to represent this elliptic solution. 



2 Movable singularities of CGL3 and CGL5 

Our interest is to count the number of distinct Laurent series for M and ip. 
The results, obtained in [3] for CGL3 and in for CGL5, are the following. 

A first set of poles xi = £ — £i — > for both M and ip arises by balancing 
A xx and the highest nonlinearity (\A\ 2 A for CGL3, \A\ 4 A for CGL5), 

CGL3 : A ~ A oX i 1+%a , A ~ A aX i^ ia , (1 - ta)(2 - ^a)p + A^g = 0, (13) 

CGL5: A~^ oX r i+iQ , A ^ A oX ^^ a , {\ - ia){^ - ia)p + A 4 r = 0, (14) 

they define two real values a and two real values Aq (CGL3) and Aq (CGL5), 

CGL3 : d t a 2 - 3d r - 2d t = 0, A\ = 3a/ di, (15) 
CGL5 : 4e,a 2 - 8e r - 3e 4 = 0, Aq = 2a/e l: (16) 

and the Fuchs indices are —1,0 and two irrational numbers 



CGL3 : indices = - 1, 0, (7 ± y/l - 24a 2 )/2, (17) 
CGL5 : indices = - 1, 0, (5 ± a/I - 32a 2 )/2. (18) 

For CGL3, M presents double poles and ip simple poles, 



M = moXi 

. dilTlo _! 

i> = — 5— xi 



1 + j Xi + 0(xi)\ , m = Al (19) 

l + £xi + 0(X?)l. (20) 
6 J 
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and for CGL5 M and tp present simple poles, 
i 2rf r r?io — 2eidirriQ 



M = moXi 1 



1 



eim —1 



4(l + e frn4) 
4(ij + 5ei<i r 



xi + o(x?) 



s 



m - 



4(l + ^mg) 



, m = Aq, (21) 
(22) 



in which both invariances (|ll[) - (112l) require changing mo to —rriQ. The number 
of distinct Laurent series M near xi is two (CGL3) or four (CGL5), and the 
number of series ip is two (CGL3), four (CGL5 q ^ 0) or two (CGL5 q — 0). 
A second set of singularities is easier to compute from the system © , 



1 

M 

if> = 



Mi 





3 -1 
2%2 



Xi 



1 + M lX 2 + < M{ 



KsMx - -g r 



7.9' 



X2 



1 + (m + Mi) X 2 + <yM( + 2«iMi + -ft 
3j«? (3di-irf r )M , 



+ o 1 (9r+j9i)ni + 



4j 

(11K 



5 . 

f 6^ 
f 4g r 



(23) 



X2 



4jft)Mi 



4 4 4 

+ 3 JKi Mi 2 + jM? } xl + G(xt)] , (24) 

in which Mo, Mi are arbitrary constants, and = —1. Invariances (|lip - (fT2"]l 
require changing Mo to —Mo, with Mi = when «j = 0. This defines either 2N 
(when q 0) or N (when q = 0) simple poles of ip, with N an undetermined 
integer. A direct study [5] of the third order ODEs for M and ip shows that 
neither M nor tp admit other movable poles. 

For CGL3 (resp. CGL5), M and ipM admit two (resp. four) Laurent series. 



3 Results from Clunie's lemma 

For convenience, £ will be denoted as z in this section only. We shall prove 

Theorem 1 For all values of the constants p,q,r,j,c,UJ, all meromorphic 
traveling wave solutions M of CGL3 and CGL5 equations belong to class W . 

The method we use is a refinement of Eremenko's method developed in 
[TU] as well as [TTirnZl IT] , based on the local singularity analysis of the solutions 
of the given differential equation and on the zero distribution and growth rate 
of their meromorphic solutions by using Nevanlinna theory. 

Several partial results have been previously obtained [2"4"1[T71I25] for finding- 
solutions of ([8]), but they are incomplete and Theorem [T] settles the question. 
These previous results are the following. 

1. For CGL3, when d r ^ 0, all solutions belonging to class W have been 
found [24j : there are six distinct solutions which are rational functions in 
one exponential function and there is no elliptic solution. 

2. For CGL3, when K; 7^ 0, there exists no elliptic solution [T7] . 

3. For CGL5, when k; = 0, there exists exactly one elliptic solution :2!) . 
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Let us recall a definition. For a differential polynomial of /, 

where j = (jo, ■■■,jk) is a multi-index and / and a,j are meromorphic functions, 
the sum jo + ---+jk is called the degree of the monomial djp° (f'Y 1 ■ ■ ■ (f^Y k ■ 
The total degree of P(z, f) is defined as the maximum of the degrees of its 
monomials. 

We shall assume the readers are familiar with the terminology and results 
of Nevanlinna theory |15, 21,25 (see [11] for a quick introduction). Here, we 
recall some basic notations of Nevanlinna theory. Let / be a non-constant 
meromorphic function on the open disc D(r) = {z : \z\ < r} where r. De- 
note the number of poles of / on the closed disc D(r) by n(r, /), counting 
multiplicity. Define the integrated counting function N(r, f) by N(r, f) = 
n(0, /)logr + JT [n(t, /) — n(0, /)] * and the proximity function m(r,f) by 

m(r, f) — / Q 2lr log + f(re t9 )^, where log + x = max {0, log x). Finally, the 
Nevanlinna characteristic function T(r, /) is defined by T(r, /) = m(r, /) + 
N(r, f) and we let S(r, f ) be a term such that x(rf) ~ * ® > as r ~ * +°°- 
To prove Theorem [TJ we make use of the well known Clunie's Lemma, 

Lemma 1 f21l 2.4-2] Let f be a transcendental meromorphic solution of 

f n P(z,f) = Q(z,f) } (25) 

with n a positive integer, P and Q differential polynomials of f with mero- 
morphic coefficients a\ such that m(r,a\) = S(r,f). If the total degree of Q 
is less than or equal to n, then 

m(r,P(z,f))=S(r,f). (26) 

Actually, all we need is the following corollary of Clunie's Lemma. 
Corollary 1 Let f be a transcendental meromorphic solution of the ODE 

r +1 =Q(z,f), (27) 

with n a positive integer and Q a differential polynomial of f with meromorphic 
coefficients a\ such that m(r, a\) — S(r, f). If the total degree of Q is less than 
or equal to n, then f must have infinitely many poles. 



Proof of Corollary [TJ 

Taking P(z,f) — f in Lemma [TJ we conclude that m(r,f) — S(r,f), and 
therefore (1 — o(l))T(r, /) = N(r, f). Assume that / has finitely many poles. 
Then N(r, f) — O(logr), and therefore T(r, /) = O(logr), which is impossible 
since / is transcendental. □ 
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Proof of Theorem [T] 

Let M be a solution of ((8]) which is meromorphic in the complex plane. If 
M is rational, then we are done. So suppose M is transcendental and let us 
prove that M has infinitely many poles. We first rewrite the second equation 
of © as 



It follows easily that if ip has infinitely many poles, then so does M. 

We first show that if ip is transcendental, then ip has infinitely many poles 
and hence so does M. One can build an ODE for ip(z) via the elimination of 
M between the system (J6j) . This third order ODE is given as follows, 



where the differential polynomial Q(z, ip) has the total degree 19 (CGL5 case) 
or 9 (CGL3 case). Applying Corollary [T] to ([29} and (j3"0"j) . we conclude that ip 

and therefore M must have infinitely many poles. 

Now suppose ip is rational, then it is well known that T(r,ip) — O(logr) 
and T(r,ip') — O(logr), and therefore m(r,ip) — S(r,f) and so is m(r,ip') 
since M is transcendental. Now (|28|) can be written as e^M 3 = —ip'M—ipM + 
(Kitp)M — diM 2 + g r M and applying Corollary [1] to it, we again conclude that 
M has infinitely many poles. 

Secondly, knowing that the transcendental meromorphic solution M has 
infinitely many poles, let us prove that it is a periodic function. By the local 
singularity analysis (section |2J), if Zo is a pole of M, CGL3 (resp. CGL5) 
admits exactly two (resp. four) Laurent series M with poles at z — zq obeying 
the ODE ©. Now let z j} j = 1,2,3, •• • be the poles of M(z); the functions 
Wj(z) = M(z + Zj — zq) are then meromorphic solutions of the ODE (jSJ) with 
a pole at zo, therefore some of them must be equal. Consequently, M is a 
periodic function. 

Without loss of generality, we assume that M has a period of 2iri. Let 
D = {z : < Imz < 2tt}. If M has more than three (CGL3) or five (CGL5) 
poles in D, then by the previous argument we conclude that M is periodic in 
D and therefore is indeed an elliptic function and we are done. 

Now suppose M has at most two (CGL3) or four (CGL5) poles in D. 
Since M is a periodic function with period 2iri, we have N(r,M) = 0(r), as 
r -> oo. It follows from (1 - o(l))T(r,M) = N(r, M) that T(r,M) = 0(r). 
By Nevanlinna's First Fundamental Theorem, we know that for any a G C, 
N(r, 1/(M — a)) = 0(r) as r — > oo. By the periodicity of M, we conclude that 
M takes each a finitely many times in D. Hence, the function R(z) — M(ln z) is 
a single- valued analytic function in the punctured plane {z:0<|z|<oo} and 
takes each a finitely many times. It follows that is a removable singularity of 
R, and R must then be a rational function. Therefore, M(z) = R(e z ) belongs 
to class W. □ 



(tpM)' - Ki(ipM) + e;M 3 + d.M 2 - g r M = 0. 



(28) 



CGL5 : e 2 (5e 2 + e 2 ) 2 ^ 20 
CGL3 : d^d* + d 2 r )iP w = 



: Q(z,ip), 



(29) 
(30) 
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4 A method to determine all solutions in class W 

Consider an iV-th order autonomous algebraic ODE, 

E(u iN) ,...,u',u) = 0, ' = d/dx, (31) 
admitting at least one Laurent series 

+oo 

u = x p ^2 UjX 1 , X = x-x . (32) 

3=0 

There exists an algorithm [24] to find in closed form all its elliptic or 
degenerate elliptic solutions. Its successive steps are [3[S]: 

1. Find the structure of movable singularities (e.g., 4 families of simple poles). 
For each subset of families (e.g. 2 families of simple poles) deduce the 
elliptic orders m, n (e.g. m = 2, n = 4) of u, v! and perform the next steps. 

2. Compute slightly more than (m + l) 2 terms in the Laurent series. 

3. Define the first order m-th degree subequation 

m 2m — 2k 

F(u, «') = 2J X! a hkU j u' = 0, a , m ^ 0. (33) 

According to results of Briot-Bouquet and Painleve [B] , any solution of (|3"Tj) 
in class W must obey such an ODE, called a "subequation" because (next 
step) it admits (13"TT) as a differential consequence. 

4. Require each Laurent series (|3"2"]l to obey F(u,u') = 0, 

F= V" !> " (E F ^ | . V; : F, =0. C-!4) 

and solve this linear overdetermined system for a^fc. 

5. Integrate each resulting ODE F(u, u') = 0. 

A similar method has later been developed [9], which also takes advan- 
tage of the Laurent series and directly searches for a canonical closed form 
representation of the elliptic solutions and their degeneracies. 

Theorem [1] implies that the subequation method is indeed able to find all 
the meromorphic traveling wave solutions of the CGL3 and CGL5 equations. 

For CGL5, the subequation method has produced a new elliptic solution 
[8], characterized by the first order, fourth degree, genus one ODE 

72 2 4 3*e^ 
F 4 = M - 2 Ki MM' 3 + — ei M' 2 (eiM 2 - 12e ) + ^ 

e* ef 

648e 2 1 ^ 

+ — j 1 (288e^ + 24 ei e M 2 - e 2 M 4 ) - —g — M 2 (aM 2 - 48e ) = 0, (35) 

3 

k 2 = 48ei, g r = 36e , e r = d r = di = 0, g % = -— ( 36 ) 

Id 
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5 Integration of subequation (|35p 

Let us first recall the differential equations of Weierstrass 

p' 2 = 4p 3 - g 2 p - 53 , C = -P, (log a)' = (. (37) 

Apart the representation as a rational function of p and p', 

polynomial(p) + polynomial(p)p / 
polynomial(p) 

elliptic functions have two main decompositions, either as a sum 
n / m \ N 

c + E h c(e - ^) + E W fc) (s - . E = °< ( 39 ) 

J=l V fc=o / j=l 

in which C,rj,a,j,Cj t k are complex constants (oj distinct), or as a quotient of 
two products of an equal number of entire functions a, 

p p 
constant f[ ^'"r] ' E a i " & = °' ( 40 ) 

in which ay, /3y are not necessarily distinct complex constants. 

To obtain the complex amplitude A, which is not elliptic, one can either 
compute the couple (M, ip) then perform the quadrature J ip d£, or compute 
the logarithmic derivative B, Eq. (|10[) . then perform the quadrature J B d£. 

By elimination with (J7J, one first deduces the real subequation for ip, 

Kiip' 4 — 4niip' 3 (Kiip + 24eo) 

+8i// 2 ( - Ki(27e? - 324e 2 ) + 1440eie ^ + 27K i ei?/' 2 + 16e ifj 3 + ^V> 4 ) 
/ 1 32 

+161 - -KiV> 8 - — eoV' 7 - 26K i eiV' 6 - 1632e 1 e V' 5 - (477e 2 + 552e^) ^ 4 
v o o 

-288 (I65e 2 + 4eg) ip 3 + m (2106e 2 - 31320c 2 ,) ip 2 

+ 2 7 3 6 (e 2 - Ael) ei e ip + 243 (-9e 4 + 56e 2 eg - 144c 4 ,) ) = 0, (41) 
then the complex subequation for i? as defined Eq. (fTU)l . 
(25' + HiB + 24ie ) (-B' - KjS - 24ie ) 2 

+2~ n (16(4B 3 - 3 Ki B 2 ) - 9(k 2 + 64ie )(4B + k;)) 2 = 0. (42) 

The degree of subequation (|4"Tj) drops from four to two when m = 0. As to 
(|4"2l . it has degree three and therefore belongs to the so-called trinomial type 
integrated by Briot and Bouquet [2j §250-251 p. 395]. 

Let us derive decompositions (|39l) or (|4U| for the solution of genus one 
equations (|3"5j) for M, (l4"Tj) for or (|42l) for £?. Three steps are required. 
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The first step is to represent the solution as a rational function of p(£ — £o) 
and p'(£ — £o), in which £o is arbitrary, and to write it in the canonical 
form (|38|) . Because of the existence of an addition formula for p, 



U , N , \ / \ 1 f P'( x l) ~ P'( x 2) 

Vxi,x 2 : p{xi + x 2 j + p(xi) + p{x 2 ) = - — — _ — - 



,(43) 



such a canonical form (1381) is not unique, and general algorithms may yield 
messy expressions by performing a shift on £o- For instance, with the (oth- 
erwise powerful) command Weierstrassf orm [16 of the computer algebra 
language Maple [22], which applies to any genus one equation, the ODE 



v! = vT - u A + u z +u + 7 
is integrated as the second degree rational function 



12p + 24VV 



22 



' 144p 2 - 24p - 251 ' 92 3 ' 93 



277 

432' 



(44) 



(45) 



while a first degree rational function cq + ci/(p — C2) is sufficient. The 
same occurs with the algorithm of Briot and Bouquet [2J §250-251 p. 395] 
to integrate binomial or trinomial equations: with (|42p . instead of yielding 
a second degree rational function rational in the fixed constants (see ([50)1 
below), it yields a third degree rational function algebraic in the fixed 
constants. Consequently, the practical method used here is to determine 
the smallest degrees of the three polynomials in ([38l) , then their coefficients 
by identification. One thus finds for the solution M of (|3"5j) . 



M = 



P 2 M = P< 



8iV (3ei + 4je )(p - ei )[3eip 2 + 4(3e 2 + 4e 2 )p + 4ei(3e 2 + 5e 2 )] 



24(3e 



No 



P ■ 

p" 

92 



4je )(p 
4(ei + je )p + 
324je? 



- ei)(p 2 - 2eip 
4e 2 — 4jeie + 



-8e 2 - 
12e 2 , j 2 



12e 2 ) + Ki P 2 M p 
= -1, 



e,(3ei +4je )' 
p(£- £0^2,53), 
= 4(p-ei)(p 2 + eip 
-24(e 2 + 2e 2 ), ff3 = 



(46) 



f 7e 2 - 
4(7e 2 



12^), 
12e 2 ) ei . 



This expression will simplify greatly as (jBT))) . Because of the correspondence 
([7]) , the solution -0 of (|4Tj) involves the same square root j of — 1 as in (|46|) . 
When Ki is nonzero this is 



P. 



i> = - 

p>i> _ 

- 945e? - 1434je?e - 1192e^e 2 - 1440jegei - 
Q2 = 9jei(ei(p 2 + 22p ei + 24jpe ) + 121e? + 48eie^ 
p:-p(e-Co M ,G 2 ,G 3 ), 

, ' 2 
P 



Qip' 



j/Cj(9ei - 4je ) 

24ei 12e x (3eip + 15e 2 + 16e 2 )((p + 2 ei ) 2 + 3(3ei + 4je ) 2 ) ' 

-j«i(3ei + 4je )[(3ei + 2je )((9e 1 - 4je )p 2 + 2(-9ei - 44je )eip) 



384e 



OJ: 



p(£ 

= 4(p + 2ei)(p 2 - 2eip - 35e 2 - 48e 2 ), 



192j e feo + 128jeg), 



(47) 



G 2 = 12(13ef + 16e5), G 3 = 8(35e 2 + 48eS)ei, 
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while for Kj = it is 

Ki = : V = \AijV3eo f 1 + „ , (48) 

or simply (but with yet another (72), 

«i = 0: ^=y^(f-^,-768e§,0). (49) 

Finally, the solution B of (|l2"jl is expressed as, 

_ «i 6Kj(3ei + 4ze ) 2 + (p + 2 ei + 3(3ei + Ue ))p' 
~ 2 2((p + 2 ei ) 2 + 3(3 ei +4je ) 2 ) ' (50) 

p:=p(f-^»G a) G 8 ), 

The properties of the above three expressions are: all coefficients (except 
the global factor No) are rational in (Ki,g r ), the two different Weierstrass 
functions p(.,<?2><?3) and p(., G2,Gs) are linked by a Landen transforma- 
tion (Appendix A), the relation (fTU)) between i?, M, tp holds true when the 
square root j of —1 is equal to +i and the constant origins £o" f ,£(f are 
equal. 

The degeneracy A — implies 32 = 33 = 0, i.e. it directly defines the 
reducible subequation (3M') 4 — e|M 8 = 0, whose solutions are rational. 
Because of this, even a four-family extension of the method used in |23j 
would fail for CGL5. 

The second step is to compute the partial fraction decomposition of the 
rational functions (|38l) of the variable p, considering for a moment p' as 
a parameter. The rational function (|46)l . once converted to the canonical 
form (|38[) . admits four poles for each choice of j, and we characterize their 
affixes £M, by choosing the signs of p'(£ j 1 ^,, g 2 , gz) as follows, 

' p(#*,flii,Sa) = (-3 + 3(/ + VSj^p + (-1) V) e x /6, 
< p'(^,32,33) = (9/ + 3((-l) 1+fc -jV3)p + j 2 -V)ei^/36, (51) 
. eip 2 = 3jV3(3 ei + 4je ), j = ±i, k = 1, 2, 3, 4. 

The rational function (|47[) admits one real pole and, for each j = ±z, two 
complex poles similarly characterized as follows, 

/ rj . / P«f , G 2 , G 3 ) = -5ex - 16e 2 /(3 ei ), 
^ • 1 p'(^, G 2 , G 3 ) = -2j Ki eo(9 e 2 + 16e 2 /(9e 2 ), j = ±i, [ ' 

f p($ fc , G 2 , G 3 ) = -2ex + (-l) fe jV3(3 ei + 4je ), j = ±i, fc = 0, 1, 

\ P%t*> G2 ' G 3) = (3 - (-l) fe jV3) Ki (3 ei + 4jeo)/2. 1 J 

Finally, the two poles of (|5U| are just p[Qu, G2, G3). 
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Modulo the periods of p(., G2, G3), the affixes of these poles obey 

Co + & - C : P = "2ei, p' = 0, C = Hi (half-period), (54) 

Co + ^1 : p = 7ei ' p ' = — 6j'«ie , (55) 

Co - Ci : P = ~ 5ei ' &>' = - 2 ^3 K i e o, (56) 

C* + C,u- fc - C : P = 00, p' = oo (period). (57) 

The Landen transformation maps p(^j 2,92,93) and p(C4,<?2) <7a) to 
p(Co>G 2 ,G 3 ), and maps 92, 9s) and p(Cj^, 52,^3) to p(Ci,G 2 ,G 3 ). 

Expressions (|4"6) , (|4"?1) . ([5Uf thus evaluate to the sum 

E' constant + constant p'(£ — £q) , . 
77 n 7^1 • (58) 

J=1 P(f-Co)-P(60 

The i/iird siep is, using the classical identities 

P» 



C(u + v) + C(u - v) - 2C(u) 



Vu.v: i ^ ' ' P(«)-PW (59) 

C(« + «) - C(« - v) - 2C(«) - , r , 
I PW - PW 

to convert l|58p into a finite sum of £ functions. 

The result for flUl, 



oi/4 4 

VKi : M = -= ^i*- 1 (C(C - $ h ,92,to) + 92, to)) , f = -1,(60) 



- e * k=i 



clearly displays the four simple poles. 

As to the three simple pole variables B, ip, M' /M, their decompositions 
evaluate to (we abbreviate £(.,G2,G 3 ) to ((.) ) 



,9ei - 4je 



(ae-Ci) + c(Ci) 



«i ^ : </> = -j + § (Ctt ~ C) + C(C) - C(0)) (61) 

+ ^ (c« - Co) + c(Co) - c« - Ci) - c(Ci)) > 

- (c(€ - Co) + C(Co) + C(€ - Ci) + c(Ci))) ■ 
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The choice j = +i must be made for (fTU)) to hold true, while the choice j = —i 
corresponds to the relation 

(. cs r ^\ 
Ae ^j=_- # . (62) 

Before taking the quadrature of the above three expressions ([BT]) . let us 
recall the definition of the element simple |14[ vol. II, p. 506] introduced by 
Hermite for integrating the Lame equation, 

H (£, q, k) = g||±gh (*-C(g))*, (k, q) constants. (63) 

Its only singularity is a simple pole with residue unity at the origin. 
Equations (jBTjh and (|6"TjU then integrate as 

— + « — _ _ 

A = K e 2p H(e! _^ 0)0) (- 1+ ^/2 Hfe _^ i)0) (- 1 - i V3)/2 ; (64) 

3ei + 4je 

0:M = K ie 12ei M H (£, 0) H(£, 0)" 1 H (£, -$ x , 0)" 1 . (65) 

The integration constants Kq,Ki arc determined by requiring that, near the 
simple pole Xi — £~£/o ^ of M, the variables A and M admit the principal 

parts A ~ A x[~ 1+iV * )/2 , M ~ A 2 oX i\ A s = 3/ef , see (HU and (ZD. 

In order to check that the product of the complex amplitude A |64|) by its 
complex conjugate is equal to the decomposition (|65p . one must take account 
of (|57|) and remember that the origin of £, not displayed in the above formulae, 
depends on j and is therefore different for A and its complex conjugate. 

The restriction k; =^ in (|65[) is removed by taking into account the relation 

C($o) + C(&) - C(#) = C($o + $,i " + «i/ 4 + J«ie /(3 ei ), (66) 
and using the definition (|5~4"|) . yielding 

„ . _,,.-*{ ^-<rwo (67) 

In order to check the equality of the two decompositions of M as the sum 
([BH)) and the product ([BT]) . it is sufficient to convert the elliptic function (|6"T|) 
to a rational function of p(£, G*2, G3) and its derivative, then to identify it to 
(|46[) modulo the Landen transformation (Appendix A). 

Numerical simulations with periodic boundary conditions [27l Fig. 4] do 
display solutions M having a real period (similar features are observed in 
CGL3 [U Fig. 7]), these could well correspond to the present elliptic solution. 

Remark. The elliptic (hence singlevalued) nature of d\og(Ae tu ' t )/d^ ex- 
plains the so-called "ad hoc Hirota method" in which A is essentially 
assumed to be a product of powers of entire functions, the powers being those 
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of the singularity structure, here (—1 ± iy/3)/2. In order to recover our result 
(|M]h two upgrades to this method are needed: (i) to assume A to be a product 
of powers of Hermite's simple elements f|63[) . not of Weierstrass a functions or 
Jacobi 9 functions, so as to ensure that the logarithmic derivative of A is ellip- 
tic; (ii) to allow arbitrary shifts £j in the arguments of the entire functions, not 
only half periods like with the choice 0j(£),j — 0, 1, 2, 3 in the Jacobi notation. 



Appendix A. Landen transformation 

We are indebted to the grateful indications of Yuri Brezhnev for this appendix. 

The Landen or Gauss transformation consists in halving only one of the two periods, it 
is naturally defined |18 | |19l p. 384] in the notation of p displaying the two periods 2u), 2a/, 

p(x\w, 2u/) = p(x\2u), 2d 1 ) + p(x - w|2w, 2ui') - p(o;|2uj, 2w'). (68) 

In the other usual notation 



(a;,G 2 ,C?3) = p(ir|oj,2a/), pi{x, g 2 , 93) = p(x\2iv, 2uj'), 
p' 2 = 4(p 3 - g 2 p - g 3 ) = 4(p - ei)(p - e 2 )(p - e 3 ), (69) 
H = 4(p? - G2P1 - G 3 ) = 4(pi - £i)(pi - E 2 )(pi - E 3 ), 



the expression of p(x\uj,2uj') as a rational function of p{x\2ui,2ui') is 

p(x,G 2 ,G 3 ) = p(x,g 2 ,g 3 ) - - 92 ~ ^ -, (70) 

4(p(a;,32,ff3) - ei) 

and similarly at the £ and <r levels 1 191 Eqs. (16b), (17b)] 

C(x,G 2 ,G 3 ) = C(z,92,33) + C(x - u,g 2 ,g 3 ) - eix + C(",92,S3), (71) 

«T(*,G a ,G a ) = e ei^ 2 /2 - C(^,<?2, 93)^ ^,32,33)^+^,32,93) ^ (72) 

o"(^,92,33) 

Between ej,Ej (and gk,Gk), there exist two algebraic relations 
£1 = -2ei, (B 2 - E 3 f = 36e 2 - 4(e 2 - e 3 ) 2 , 

-329233 + 2293G2 + H92G3 - G2G3 = 0, (73) 
196s! + 49g^G 2 - 7260c; 2 + 660g 3 G 3 - 15G§ = 0. 

The ratio —2 of the two zeros — 2ei of p' x 2 in II47I I and e\ of p' 2 in (146 I t is the signature of 
such a Landen transformation. 
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